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A METHOD AND AN APPARATUS FOR PREDICTING INTAKE 
MANIFOLD PRESSURE OF AN INTERNAL- COMBUSTION ENGINE 

BACKGROUND OF THE INVENTION: 
Field of the Invention: 

[0001] The present invention relates to a method and an apparatus for 
predicting an intake manifold pressure of an internal-combustion engine . In 
particular, the present invention relates to a method and a fuzzy estimator 
for predicting an intake manifold pressure, using fuzzy algorithm. 

Description of the Related Art: 

[0002] Fuel injection control is carried out for efficient combustion in 
internal-combustion engines. Fig. 1 shows an intake portion of an, 
internal-combustion engine, to which the present invention is applied. A 
drawn air is supplied to a cylinder through a throttle 1. Throttle opening is 
controlled at a desired value. An intake manifold pressure is measured with 
a sensor 2. For appropriate fuel injection control, an amount of air to be 
drawn into the cylinder must be estimated and therefore a value of intake 
manifold pressure must be predicted. 

[0003] Conventionally, such an intake manifold pressure predicting 
apparatus as described below was disclosed, for example in Publication of 
Japanese Unexamined Patent Publication (KOKAI) No. 2-42160. The 
apparatus predicts an intake manifold pressure at the present time based on 
operating states of the internal-combustion engine and predicts intake 
manifold pressure at a time of prediction, a certain period ahead, based on 
the predicted intake manifold pressure and a detected value of intake 
manifold pressure. 
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[0004] However, the above-mentioned predicting apparatus was not able 
to make a prediction with sufficient accuracy over a wide range of number of 
revolutions of the internal-combustion engine. The reason is that the 
apparatus corrects a predicted intake manifold pressure during a period 
while an intake valve is transiently closed, evenly based on an error in intake 
manifold pressure AP and independently of operating states of the 
internal-combustion engine. 

[0005] Japanese Patent No. 2886771 discloses an intake manifold 
pressure predicting apparatus which makes a prediction considering 
operating states of the internal-combustion engine to solve the 
above-mentioned problem, and enables a high accuracy control even in the 
case of lower number of revolutions and higher intake manifold pressures. 

[0006] In the conventional method mentioned above, a predicted value of 
intake manifold pressure (hereinafter referred to as HATPB) is calculated 
based on a difference (hereinafter referred to as APB) between successive 
values of intake manifold pressure (hereinafter referred to as PB) and a 
difference (hereinafter referred to as ATH) between successive values of 
throttle opening. Further, HATPB is used for fuel injection control and 
retrieval of parameters for fuel adhesion correction. Now, ATH and APB 
are represented as below provided that k is a point in control time 
synchronized with intake stroke (TDC). 

ATH(k) = TH(k) + TH(k- 1) (l) 

APB(k) = PB(k) + PB(k-l) (2) 
[0007] On the other hand, between the detecting portion of a PB sensor 
and an intake manifold, on the intake manifold side or in the PB sensor, a 
labyrinth mechanism or the like has recently been provided to prevent water 
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from entering there. Accordingly, a time lag and a time delay between an 
actual value of pressure and an output of the PB sensor, have become larger. 

[0008] There have been attempts to use conventional prediction 
algorithms to compensate the time delay. However, the algorithms have 
caused an overshoot of HATPB against an actual PB, as shown in Fig. 2 or a 
discontinuous behavior of HATPB. The reason for the overshoot is that the 
conventional prediction algorithms make up for insufficient accuracy of a 
predicted value of PB through feedback of an error between a predicted value 
of PB a certain time in advance and the current value of PB. Further, the 
reason for the discontinuous behavior is that the conventional prediction 
algorithms use one of predicted values calculated respectively based on APB 
and ATH, by switching according to certain conditions. Such behaviors of 
the conventional prediction algorithms have such an influence on fuel 
injection control as to cause a problem that variation in variables during 
transient operations becomes larger to increase an amount of emissions of 
offensive exhaust elements. 

[0009] Accordingly, there has been a need for a new prediction algorithm 
for PB to compensate for a larger lag and a larger time delay without 
producing an overshoot or a discontinuous behavior of HATPB. 

SUMMARY OF THE INVENTION: 

[0010] A method for predicting intake manifold pressure according to an 
aspect of the present invention, comprises the step of obtaining a difference of 
values of intake manifold pressure and a difference of values of throttle 
opening. The method further comprises the step of obtaining a predicted 
difference of values of intake manifold pressure, through algorithm of 
estimation with fuzzy reasoning. The method further comprises the step of 



adding the predicted difference of values of intake manifold pressure, to a 
value of intake manifold pressure, to obtain a predicted value of intake 
manifold pressure. The algorithm of estimation with fuzzy reasoning, 
includes fuzzy rules determined based on an amount of a difference of values 
of intake manifold pressure and an amount of a difference of values of 
throttle opening. 

[0011] An apparatus for predicting intake manifold pressure, according to 
another aspect of the present invention, comprises a device for obtaining a 
difference of values of intake manifold pressure, a device for obtaining a 
difference of values of throttle opening and a fuzzy estimator. The fuzzy 
estimator receives as inputs the difference of values of intake manifold 
pressure and the difference of values of throttle opening and obtains and 
outputs a predicted difference of values of intake manifold pressure, through 
algorithm of estimation with fuzzy reasoning. The algorithm of estimation 
with fuzzy reasoning includes fuzzy rules determined based on an amount of 
a difference of values of intake manifold pressure and an amount of a 
difference of values of throttle opening. 

[0012] A computer-readable medium, according to another aspect of the 
present invention, has a program stored therein, which is made to perform 
the step of obtaining a difference of values of intake manifold pressure and a 
difference of values of throttle opening. The program is further made to 
perform the step of obtaining a predicted difference of values of intake 
manifold pressure, through algorithm of estimation with fuzzy reasoning. 
The program is further made to perform the step of adding the predicted 
difference of values of intake manifold pressure, to a value of intake manifold 
pressure, to obtain a predicted value of intake manifold pressure. The 
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algorithm of estimation with fuzzy reasoning, includes fuzzy rules 
determined based on an amount of a difference of values of intake manifold 
pressure and an amount of a difference of values of throttle opening. 

[0013] Thus, in the aspects of the present invention mentioned above, a 
predicted value of an output of the PB sensor is calculated through algorithm 
of estimation with fuzzy reasoning, based an on output of the PB sensor and 
TH opening. Then, an amount of fuel injection into the internal combustion 
engine is determined based on the predicted value. In particular, use of 
fuzzy rules based on an amount of APB and that of ATH allows a control 
effectively containing information on a change in TH which is ahead of a 
change in PB. According to the aspects of the present invention mentioned 
above, even when there occurs a large time delay or a large lag between an 
actual value of intake manifold pressure (negative) and an output of the PB 
sensor, predicted values will not be discontinuous as in conventional methods 
and predicted values which are continuous can be calculated with higher 
accuracy. Accordingly, an air-fuel ratio of the internal combustion engine 
will not show a discontinuous behavior. Further, an overshoot of the 
predicted value against the actual value of intake manifold pressure, can be 
dramatically reduced compared with conventional methods. 

[0014] According to an embodiment of the present invention, a difference 
of values of intake manifold pressure is classified based on its amount into 
positive one, that of zero or negative one and a difference of values of throttle 
opening is classified based on its amount into positive one, that of zero or 
negative one. Further, fuzzy rules are provided respectively for areas 
determined by the two kinds of classifications. 
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[0015] Accordingly, operations for estimation through fuzzy reasoning 
will not become complicated and can be carried out simply and efficiently. 

[0016] According to another embodiment of the present invention, a 
second order difference of values of intake manifold pressure, is further 
obtained. Then, fuzzy rules are determined based on an amount of a 
difference of values of intake manifold pressure, an amount of a difference of 
values of throttle opening and an amount of a second order difference of 
values of intake manifold pressure. 

[0017] Accordingly, a control effectively using information on behaviors of 
APB in the future, contained in AAPB (a second order difference of values of 
intake manifold pressure), can be carried out. 

[0018] According to another embodiment of the present invention, a 
second order difference of values of intake manifold pressure is classified 
based on its amount into positive one, that of zero or negative one. Then, 
fuzzy rules are provided respectively for areas determined by three kinds of 
classifications based on an amount of a difference of values of intake manifold 
pressure, an amount of a difference of values of throttle opening and an 
amount of a second order difference of values of intake manifold pressure. 

[0019] Accordingly, operations for estimation through fuzzy reasoning 
will not become complicated and can be carried out simply and efficiently. 

[0020] According to another embodiment of the present invention, a value 
obtained by delaying a throttle opening value by a time delay, is used as a 
throttle opening value. 
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[0021] Accordingly, even if a volume of the intake manifold is so large that 
there occurs a time delay between a change in throttle opening and a change 
in actual PB, information on TH can be used for prediction and predicted 
values which do not shown discontinuous behaviors, can be calculated with 
higher accuracy. 

[0022] According to another embodiment of the present invention, a 
relationship between a throttle opening value and a desired value of throttle 
opening is modeled using a time delay element and a lag system and a value 
estimated through the model and the desired value is used as a throttle 
opening value. That is, provided that an estimated value of throttle opening 
and a desired value of throttle opening at point k in time, are respectively 
represented as THHAT(k) and THCMD(k), a value corresponding to a time 
delay is represented as ddly and a constant is represented as Kdly, an 
estimated value of throttle opening THHAT(k) is obtained by the following 
equation for use. 

THHAT(k) = Kdly x THHAT(k) + (l - Kdly) x THCMD(k - ddly) 
[0023] Accordingly, even when an electronically controlled throttle 
involving a time delay before reading of TH, is used, a predicted value can be 
calculated with high accuracy. 

[0024] According to another embodiment of the present invention, the 
membership function for the consequent part of the algorithm of estimation 
with fuzzy reasoning, is a bar-shaped singleton function. 

[0025] Accordingly, operations for estimation through fuzzy reasoning 
will not become complicated and can be carried out simply and efficiently. A 
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processor proof against such service conditions for vehicles as cryogenic 
temperatures, high temperatures, high humidity and vibration, can hardly 
be provided with such high computing ability as is able to carry out mini-max 
gravity method of the algorithm of estimation with fuzzy reasoning. 
However, use of a bar-shaped singleton function as the membership function 
for the consequent part, enables the processor to carry out mini-max gravity 
method and to calculate a predicted value with high accuracy. 

[0026] According to another embodiment of the present invention, inputs 
are subjected to filtering. 

[0027] Accordingly, use of filtered data as input data to the algorithm of 
estimation with fuzzy reasoning, prevents a predicted value from oscillating 
even when noises are mixed. 

[0028] According to another embodiment of the present invention, the 
filtering is carried out with an adaptive filter. 

[0029] Accordingly, noises can be eliminated from the data to a sufficient 
degree while maintaining a phase delay of the data minimum so that the 
data can be used for predicting operations. 

[0030] A method for obtaining a predicted value of a variable, according to 
an aspect of the present invention, comprises the step of obtaining a 
difference of values of a variable to be predicted and a difference of values of 
another variable ahead of the variable to be predicted. The method further 
comprises the step of filtering the differences with adaptive filters. The 
method further comprises the step of obtaining a predicted difference of 
values of the variable to be predicted, through algorithm of estimation with 
fuzzy reasoning. The method further comprises the step of adding the 
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predicted difference of values of the variable to be predicted, to a current 
value of the variable to be predicted, to obtain a predicted value of the 
variable to be predicted. The algorithm of estimation with fuzzy reasoning 
includes fuzzy rules determined based on an amount of a difference of values 
of the variable to be predicted and an amount of a difference of values of the 
variable ahead of the variable to be predicted. 

[0031] A predicting apparatus, according to an aspect of the present 
invention, comprises filters for filtering inputs and a fuzzy estimator. The 
fuzzy estimator receives as inputs a difference of values of a variable to be 
predicted and a difference of values of another variable ahead of the variable 
to be predicted and obtains and outputs a predicted difference of values of the 
variable to be predicted, through algorithm of estimation with fuzzy 
reasoning. The algorithm of estimation with fuzzy reasoning includes 
fuzzy rules determined based on an amount of a difference of values of the 
variable to be predicted and an amount of a difference of values of the 
variable ahead of the variable to be predicted. 

[0032] A computer-readable medium, according to an aspect of the 
present invention, has a program stored therein, which is made to perform 
the step of obtaining a difference of values of the variable to be predicted and 
a difference of values of another variable ahead of the variable to be 
predicted. The program is further made to perform the step of filtering the 
differences with adaptive filters. The program is further made to perform 
the step of obtaining a predicted difference of values of the variable to be 
predicted, through algorithm of estimation with fuzzy reasoning. The 
program is further made to perform the step of adding the predicted 
difference of values of the variable to be predicted, to a current value of the 
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variable to be predicted, to obtain a predicted value of the variable to be 
predicted. The algorithm includes fuzzy rules determined based on an 
amount of a difference of values of the variable to be predicted and an amount 
of a difference of values of the variable ahead of the variable to be predicted. 

[0033] Thus, according to the aspects mentioned above, use of fuzzy rules 
based on an amount of a difference of the variable to be predicted and that of 
the variable ahead of the variable to be predicted, allows a control effectively 
containing information on a change in the variable ahead of the variable to be 
predicted. Further, use of adaptive filters allows noises to be eliminated 
from the data to a sufficient degree while maintaining a phase delay of the 
data minimum so that the data can be used for predicting operations. 

BRIEF DESCRIPTION OF THE DRAWINGS: 

[0034] Fig. 1 shows an intake portion of an internal-combustion engine. 

[0035] Fig. 2 shows behaviors of a value predicted by conventional 
algorithms for predicting intake manifold pressure. 

[0036] Fig. 3 shows relationships among PB, TH, APB, ATH and AFZPB. 

[0037] Fig. 4 shows membership functions for the antecedent part 
according to an embodiment of the present invention. 

[0038] Fig. 5 shows membership function for the consequent part 
according to an embodiment of the present invention. 

[0039] Fig. 6 shows fuzzy rules used in an embodiment of the present 
invention. 

[0040] Fig. 7 shows a state to which Rule 1 is applied. 
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[0041] Fig. 8 shows a state to which Rule 2 is applied. 

[0042] Fig. 9 shows a state to which Rule 3 is applied. 

[0043] Fig. 10 shows a state to which Rule 4 is applied. 

[0044] Fig. 11 shows a state to which Rule 5 is applied. 

[0045] Fig. 12 shows a state to which Rule 6 is applied. 

[0046] Fig. 13 shows a state to which Rule 7 is applied. 

[0047] Fig. 14 shows a state to which Rule 8 is applied. 

[0048] Fig. 15 shows a state to which Rule 9 is applied. 

[0049] Fig. 16 shows a method by which the degree of fulfillment for Rule 
6 is obtained based on the membership functions for the antecedent part. 

[0050] Fig. 17 shows a method by which a weight of Rule 6 is obtained 
based on the membership function for the consequent part. 

[0051] Fig. 18 shows the result of estimation of FZPB using the fussy 
inference algorithm according to an embodiment of the present invention. 

[0052] Fig. 19 shows the membership functions of the antecedent part 
according to another embodiment of the present invention. 

[0053] Fig. 20 shows fuzzy rules according to another embodiment of the 
present invention. 

[0054] Fig. 21 shows the membership function of the consequent part 
according to another embodiment of the present invention. 
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[0055] Fig. 22 shows a state to which Rule 10 is applied. 

[0056] Fig. 23 shows a state to which Rule 11 is applied. 

[0057] Fig. 24 shows a state to which Rule 6 is applied according to 
another embodiment of the present invention. 

[0058] Fig. 25 shows a state to which Rule 9 is applied according to 
another embodiment of the present invention. 

[0059] Fig. 26 shows the result of estimation of FZPB using the fussy 
inference algorithm according to another embodiment of the present 
invention. 

[0060] Fig. 27 shows a fuzzy estimator according to an embodiment of the 
present invention. 

[0061] Fig. 28 shows a fuzzy estimator provided with adaptive filters 
according to an embodiment of the present invention. 

[0062] Fig. 29 shows behaviors of FZPB with and without the use of 
adaptive filters. 

[0063] Fig. 30 shows a state in which an electronically controlled throttle 
is used. 

[0064] Fig. 31 shows a fuzzy estimator for the use with an electronically 
controlled throttle. 

[0065] Fig. 32 shows the result of prediction of a fuzzy estimator for the 
use with an electronically controlled throttle. 
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[0066] Fig. 33 shows a fuzzy estimator for the use with a large-volume 
intake manifold. 

[0067] Fig. 34 shows a procedure of the embodiment of the present 
invention in which adaptive filters are used. 

[0068] Fig. 35 shows an example of an electronic control unit used in 
embodiments of the present invention. 

DETAILED DESCRIPTION OF THE PREFFERRED EMBODIMENTS: 
[0069] In the present invention, a predicted amount of change in PB 
(hereinafter referred to as AFZPB) is obtained through estimation with fuzzy 
reasoning, using fuzzy algorithm including fuzzy rules defined based on an 
amount of APBand that of ATH. Then, a predicted value (hereinafter 
referred to as FZPB) through estimation with fuzzy reasoning is calculated 
by the following equation. 

FZPB(k) = PB(k) + AFZPB(k) (3) 

[0070] That is, AFZPB(k) is added to the current PB sample value PB(k) 
to obtain FZPB(k). "k" indicates a point in control time synchronized with 
intake stroke (TDC). Fig. 3 shows relationships among PB, TH, APB, ATH 
and AFZPB. 

[0071] Fig. 6 shows fuzzy rules used in an embodiment of the present 
invention. APB is classified based on its amount into positive one, that of 
zero or negative one. ATH is classified based on its amount into positive one, 
that of zero or negative one. Fuzzy rules are provided respectively for 9 
areas determined by the two kinds of classifications. It should be noted here 
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that a change in TH is ahead of a change in PB and therefore contains 
information on behaviors of PB in the future. Accordingly, use of fuzzy rules 
based on an amount of APB and that of ATH allows a control effectively 
containing information on a change in TH which is ahead of a change in PB. 

[0072] Figs. 4 and 5 respectively show membership functions for the 
antecedent part and that for the consequent part. The membership 
functions for the antecedent part for APB and ATH are set trapezoidal for 
positive (P) and negative (N) and set triangular for zero (Z). As the 
membership function for the consequent part, a bar-shape singleton function 
is used for simple operations of estimation with fuzzy reasoning. 

[0073] Now each rule (Rule 1 to Rule 9 in Fig. 6) in each area will be 
described with reference to drawings. 

[0074] Fig. 7 shows a state in which both APB and ATH are negative. 
Rule 1 is applied to the state. Since both APB and ATH are negative, AFZPB 
of the membership function for the consequent part is also set negative. 

[0075] Fig. 8 shows a state in which APB is negative and ATH is zero. 
Rule 2 is applied to the state. Since APB is negative while ATH which is 
ahead of APB is zero, AFZPB of the membership function for the consequent 
part is set negative. 

[0076] Fig. 9 shows a state in which APB is negative and ATH is positive. 
Rule 3 is applied to the state. The state corresponds to the case in which 
during operation of engine brake number of revolutions of the engine 
increases more rapidly than an amount of air passing through the throttle 
increases due to increase in TH. AFZPB of the membership function for the 
consequent part is set zero. 
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[0077] Fig. 10 shows a state in which APB is zero and ATH is negative. 
Rule 4 is applied to the state. Since ATH which is ahead of APB is negative, 
AFZPB of the membership function for the consequent part is set negative. 

[0078] Fig. 11 shows a state in which both APB and ATH are zero. Rule 5 
is applied to the state. Since both APB and ATH are zero, AFZPB of the 
membership function for the consequent part is also set zero. 

[0079] Fig. 12 shows a state in which APB is zero and ATH is positive. 
Rule 6 is applied to the state. Since ATH which is ahead of APB is positive, 
AFZPB of the membership function for the consequent part is set positive 
(Pe). 

[0080] Fig. 13 shows a state in which APB is positive and ATH is negative. 
Rule 7 is applied to the state. The state corresponds to the case in which due 
to an external force number of revolutions of the engine decreases more 
significantly than it would decrease due to decrease in TH. AFZPB of the 
membership function for the consequent is set zero. 

[0081] Fig. 14 shows a state in which APB is positive and ATH is zero. 
Rule 8 is applied to the state. AFZPB of the membership function for the 
consequent is set positive (Ps). Since ATH which is ahead of APB is zero, Ps 
is set smaller than P6in Rule 6. 

[0082] Fig. 15 shows a state in which both APB and ATH are positive. 
Rule 9 is applied to the state. AFZPB of the membership function for the 
consequent is set positive (Pio). Since APB has already been positive, Piois 
set smaller than P6 in Rule 6 for the state in which APB is zero. 
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[0083] Final AFZPB is calculated based on mini-max gravity method 
using the membership functions and fuzzy rules mentioned above. The 
method will be described in detail for the case to which Rule 6 is applied, as 
an example. 

[0084] Now mini-max selection will be described. Current sample values 
of APB and ATH are represented as APB(k) and ATH(k). Degrees of 
fulfillment of these values for Rule 6 will be obtained. As shown in Figs. 6 
and 12, Rule 6 is applied to an area in which APB is zero and ATH is positive. 
Thus the fulfillment of APB for the membership function (zero) of the 
antecedent part is mAPB(6) as shown in Fig. 16(a). Further, the fulfillment 
of ATH for the membership function (positive) of the antecedent part is 
mATH(6) as shown in Fig. 16(b). In mini-max selection, the smallest degree 
of fulfillment among the degrees of fulfillment for the antecedent part, is 
selected as the degree of fulfillment m(i) for rule i. The following 
relationship is established. 

mAPB (6) < mATH (6) (4) 

[0085] Accordingly, for the degree of fulfillment for Rule 6, mAPB(6) is 
selected. 

m (6) = mAPB (6) (5) 

[0086] Further, a position of the membership function of the consequent 
part for Rule 6 is pP6 and a reference weight (length of bar) of that is wP6 as 
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shown in Fig. 17. Accordingly, weight of Rule 6 in estimation of AFZPB, that 
is weight w(6) at position pP, will be given as below. 

w(6) = m(6) x W P6 (6) 

[0087] Then, a weight w(i) is obtained for each rule i in a similar way. All 
the weights thus obtained are used for estimation and therefore this is "max" 
selection. 

[0088] Estimation of AFZPB, that is defuzzification of a fuzzy output will 
be carried out using a gravity method shown in the following equation, based 
on weights for the rules. 

9 

2 w(i) x w Pi x pPi 

AFZPB (k) = ifiJ — g (7) 

2 w(i) x w Pi 
i=i 

[0089] Thus, equation (3) is represented as below and AFZPB(k) can be 
calculated through it. 

9 

2 w(i) x wPi x pPi 

AFZPB (k) == PB (k) = (8) 

2 w(i) x w Pi 
i = i 
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[0090] Fig. 18 shows the result of estimation of FZPB using the fussy 
inference algorithm mentioned above. In Fig. 18 FZPB is closely analogous 
to the desirable predicted value, even though there exists some overshoot 
against the desirable predicted value. Thus, accuracy of estimation has 
significantly increased compared with that of the conventional algorithm 
shown in Fig. 2. 

[0091] Now another embodiment of the present invention will be 
described below. The prediction with fuzzy reasoning in the embodiment 
mentioned above, produces some overshoot as shown in Fig. 18. Such an 
overshoot tends to occur particularly under hard acceleration or snap. In 
order to eliminate such an overshoot, a change in APB, that is second order 
difference of PB (hereinafter referred to as AAPB) should be noted. AAPB is 
defined by the following equation where "k" indicates a point in control time 
synchronized with intake stroke (TDC). 

AAPB(K) = APB(k) - APB(K- 1) (9) 

[0092] Further, fuzzy rules will be defined based on an amount of APB, 
that of ATH and that of AAPB. Thus, the three kinds of classifications due to 
an amount of APB, that of ATH and that of AAPB, generate 27 areas, for each 
of which a fuzzy rule is provided. Fig. 20 shows fuzzy rules according to the 
present embodiment. Fuzzy rules shown in Fig. 20 are positioned in three 
dimensions, while those shown in Fig. 6 are positioned in two dimensions. It 
should be noted here that AAPB contains information on behaviors of APB in 
the future. Accordingly, use of fuzzy rules based on an amount of AAPB 
allows a control containing information on the information on behaviors of 
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APB in the future. Fig. 19 shows the membership functions of the 
antecedent part according to the present embodiment. 

[0093] Compared with fuzzy rules shown in Fig. 6, in Fig. 20 Rule 8 and 
Rule 9 alone differ respectively when AAPB is positive and when AAPB is 
negative. The other rules remain unchanged independently of AAPB. 

[0094] Now new rules (Rules 10 and 11 in Fig. 20) will be described with 
reference to drawings. 

[0095] Fig. 22 shows a state in which APB, ATH and AAPB are positive. 
Rule 10 is applied to the state. Since APB, ATH and AAPB are positive, 
AFZPB of the membership function for the consequent is set to the largest 
positive value (Pio). 

[0096] Fig. 23 shows a state in which APB is positive, ATH is zero and 
AAPB is negative. Rule 11 is applied to the state. AFZPB of the 
membership function for the consequent part is set zero. 

[0097] Fig. 24 shows a state in which APB is zero and ATH is positive, 
according to the present embodiment. Rule 6 is applied to the state. As 
mentioned above, Rule 10 is applied to a state in which APB, ATH and AAPB 
are positive. Accordingly, an area to which Rule 6 is applied is reduced, 
compared with the case in which AAPB is not used. 

[0098] Fig.25 shows a state in which both APB and ATH are positive, 
according to the present embodiment. Rule 9 is applied to the state. As 
mentioned above, Rule 10 is applied to a state in which APB, ATH and AAPB 
are positive. Accordingly, an area to which Rule 9 is applied is reduced, 
compared with the case in which AAPB is not used. 



19 



[0099] Fig. 21 shows the membership function of the consequent part 
according to the present embodiment. As mentioned above, Pio is the largest 
positive value. Compared with the membership function of the consequent 
part shown in Fig. 5, pP6 (position of P6> and pP9 (position of Pa) has 
approached zero. The reason is that areas to which Rule 6 and Rule 9 are 
applied, have been reduced. 

[0100] Fig. 26 shows the result of estimation of FZPB using the fussy 
inference algorithm according to the present embodiment. In Fig. 26 such 
an overshoot as shown in Fig. 18 has been eliminated. 

[0101] In the present embodiment when a time delay to be compensated is 
small, P6 and P9 having approached zero, can be set zero for the sake of 
simplifying data setting. Further, the state to which Rule 8 is applied and in 
which AAPB is positive, rarely occurs. Accordingly, this rule may be 
eliminated, or the consequent part may be set zero to invalidate its 
contributions to the prediction. 

[0102] Now embodiments in which an input is subjected to filtering, will 
be described. The present invention employs algorithm of estimation with 
fuzzy reasoning, using APB, ATH and AAPB as inputs. Accordingly, if noises 
are mixed into TH and / or PB, the values of difference and second order 
difference will be oscillating or have spikes. As a result, FZPB estimated 
through estimation with fuzzy reasoning will also be oscillating or have 
spikes at times. 

[0103] In the present embodiment, APB, ATH and AAPB are subjected to 
filtering and then input to a fuzzy estimator as shown Fig. 28. On the 
contrary, in Fig. 27 APB, ATH and AAPB are input to a fuzzy estimator. In 



20 



-« 



Fig. 27 devices for obtaining difference are represented with reference 
numerals 2711 to 2713 and a fuzzy estimator is represented with reference 
numeral 2701. In Fig. 28 devices for obtaining difference are represented 
with reference numerals 2811 to 2813, filters are represented with reference 
numerals 2821 to 2823 and a fuzzy estimator is represented with reference 
numeral 2801. 

[0104] As the filters for APB, ATH and AAPB, such an adaptive filter as 
represented in the following equations, is employed. This type of filters 
prevents FZPB from oscillating or having spikes due to noises as shown in 
Fig. 29. 

Xf (k) = X _ f (k - 1) + KP(k) • ide(k) (10) 

KP(k) = P(k - 1) / (1 + P (k-D) (11) 

ide(k) = X_f <k- 1) -X(k) (12) 

P(k + 1) = (1 / X 2 )[1 • X 2 • P(k) KX 1 + x 2 • PG0)] (13) 

X_f represents adaptive filter values for APB, ATH and AAPB, while X 
represents sample values of APB, ATH and AAPB. X x and X 2 represent 
weighting parameters. 

[0105] Now, another embodiment will be described, in which TH is 
determined in consideration of time delay, based on a desired value of 
throttle opening. Recently an electronically controlled throttle has been 
widely employed to fill the need for coordinated control with the transmission 
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for better fuel economy and for control for stabilizing steering. An 
electronically controlled throttle is often driven by a driver separate from the 
electronic control unit. The driver is connected to the electronic control unit 
via a network on the vehicle (CAN or the like). 

[0106] Accordingly, a period (10 milliseconds) of communication between 
the electronic control unit and the driver, causes a time delay between 
throttle opening command THCMD calculated by the electronic control unit 
and an actual throttle opening THACT of the electronically controlled 
throttle. That is, a time delay occurs before THCMD exerts an influence 
upon TH. Further, another time delay occurs before THACT exerts an 
influence upon TH observed by the electronic control unit via the CAN or the 
like. Fig. 30 shows a state in which an electronically controlled throttle is 
used. In Fig. 30 THCMD represents a command (reference) of throttle 
opening, THACT represents an actual throttle opening and TH represents an 
(observed) value of throttle opening, as mentioned above. 

[0107] Thus, a change in TH observed by the electronic control unit lags 
behind a change in PB and therefore ATH cannot be used to predict PB as in 
the algorithm of estimation with fuzzy reasoning mentioned above. 

[0108] Accordingly, in the present embodiment THHAT is estimated 
based on THCMD in consideration of a time delay due to communication and 
a delay in response of the electronically controlled throttle, as described 
below. A difference in THHAT, that is AD DHAT is used instead of ATH. 

THHAT(k) = Kdly x THHAT(k) + (l - Kdly) x THCMD(k - ddly) (14) 
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ATHHAT(k) = THHAT(k) - THHAT(k-l) (15) 
[0109] ddly represents a value corresponding to a time delay and Kdly 
represents a constant. 

[0110] Fig. 31 shows a system configuration including devices and Fig. 32 
shows the result of prediction, according to the present embodiment. In Fig. 
31 devices for obtaining difference are represented with reference numerals 
3111 to 3113, filters are represented with reference numerals 3121 to 3123, a 
module for carrying out operation represented by Equation (14) is 
represented with reference numeral 3131 and a fuzzy estimator is 
represented with reference numeral 3101. 

[0111] Now, another embodiment in which TH is delayed by a time delay, 
will be described. Recently a volume of the intake manifold of the engine in 
many cars has been made very large to fill the need for a larger torque at 
lower speeds. In these cases a change in an actual pressure at the intake 
manifold, lags a time delay dth behind a change in TH. Under the situation 
a change in TH is too ahead of a change in an actual pressure at the intake 
manifold to use TH for operations for prediction of an actual pressure at the 
intake manifold. In order to solve this problem, TH is delayed by a time 
delay dth before being used in operations for prediction. In Fig. 33 devices 
for obtaining difference are represented with reference numerals 3311 to 
3313, a time delay module is represented with reference numeral 3341 and a 
fuzzy estimator is represented with reference numeral 3301. 

[0112] A procedure of the embodiment of the present invention in which 
adaptive filters are used, will be described with reference to a flowchart 
shown in Fig. 34. In step S10 a PB sensor is checked for whether it is active 
or not. If the PB sensor is not active, the process proceeds with step S 80, in 
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which a substitute value is given for a sample value of PB. Further, in step 
S 90 the substitute value as the sample value is set to FZPB and the process 
is terminated. If the PB sensor is active, the process proceeds with step S 
20, in which the PB sensor is checked for whether it normally functions or 
not. If it does not normally function, the process proceeds with step S 80. If 
it normally functions, the process proceeds with step S 30, in which an 
estimated value THHAT of TH is calculated (Equation (14)). THHAT is 
used instead of TH for the case of an electronically controlled throttle. Then, 
the process proceeds with step S40, in which APB, AAPB and ATHHAT are 
calculated (Equations (2), (7) and (14)). Further, the process proceeds with 
step S50, in which operations for adaptive filters are carried out (Equations 
(10) to (13)). Then, in step S60 AFZPB is estimated. In step S70 AFZPB is 
added to a sample value of PB to obtain FZPB and the process is terminated. 

[0113] An example of an electronic control unit used in embodiments of 
the present invention, will be described with reference to Fig. 35. The 
electronic control unit includes a CPU 3501, a ROM 3511, a flash memory 
3512, a RAM 3513, an I/O unit 3514 and a communication controller 3515 for 
a network on the vehicle. The above devices are connected with one another 
via a bus 3520. 

[0114] Algorithm for predicting intake manifold pressure according to the 
present invention, may be stored as a program in the ROM 3511 or the flash 
memory 3512. Some part of the algorithm, for example fuzzy rules, may be 
stored in the flash memory 3512, while the other part may be stored in the 
ROM 3511. Alternatively, the algorithm may be stored in another type of 
memory not shown in the drawing. 



24 



